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INTRODUCTION 


A gun-launched  sitenied-range  guided  projectile  ( FR J:  ) 
is  confronted  with  a transitional  problem  in  that  it  Is 
launched  from  a rifle  barrel  as  a rotating  bullet-shaped 
member  and  must  then  convert  into  an  almost  non-rotating 
finned  missile  before  firing  a self-contained  rocket  propul- 
sion system.  Tne  resident  guidance  system  places  a con- 
straint upon  the  roll  rate  residual  of  approximately  +/-  one 
(1)  revolution  per  second  after  tne  fins  have  opened  and 
provided  roll  damping.  Since  the  residual  roll  rate  of  the 
projectile  is  an  open-loop  situation  (i.e.,  no  roll  control 
is  provided),  there  are  two  principal  difficulties  associ- 
ated with  meeting  a design  goal  on  the  residual  roll  rate. 
One  problem  is  relative  to  fin  alignment,  vith  the  residual 
roll  rate  tolerance  requiring  the  fins  to  be  aligned  with 
the  freestream  to  a small  fraction  of  a degree.  Another 
problem  concerns  the  question  whether  fin-opening  dynamics 
introduces  stresses  into  the  fin  such  that  structural  yield- 
ing occurs  thereby  providing  a separate  source  for  fin 
angular  misalignment,  and  as  a matter  of  fact,  possibly 
negating  the  time-consuming  efforts  spent  initially  in  fin 
alignment . 


The  case  considered  in  this  report  will  be  the  develop- 
ment of  fin-opening  dynamics  for  an  eignt-inch  ERJP  with  a 
six  fin  configuration.  Although  the  muzzle  velocity  of  the 
KKIP  is  fairly  well  defined  from  the  typical  43  millisecond 
acceleration  period  in  the  gun  barrel,  tne  initial  condition 
upon  free-flight  roll  rate  is  a variable  with  maximum  values 
being  on  the  order  of  33  revolutions  per  second  prior  to  fin 
opening.  The  initial,  fin  retracted  roll  rate  is  dependent 
upon  the  configuration  of  the  obturator  rings,  which  are  in- 
stalled on  the  ERDP  to  act  as  slip  rings  between  the  rifle 

barrel  and  t ft.  S miSSilc.  T u 6 ObtiirdtOT  1 1 fi  5 S g 0 t ii  f G U £ ii  a fi 

engraving  process  with  the  rifling,  and  it  is  the  relative 
slippage  between  the  rings  and  the  missile  which  prevents  a 
complete  transfer  of  angular  roll  momentum  duri  n«,  the 
initial  firing  process. 


The  ER1P  missile  may  be  assumed  as  having  a constant 
velocity  along  its  _\j.s  cf  symmetry  wnile  conserving  angular 
momentum  during  the  t lr,-op°- lag  process.  The  dissipative 
action  of  the  fins  upon  the  roll  rate,  due  to  aerodynamic 
damping  in  roll,  is  present  to  some  extent  during  fin  open- 
ing but  will  be  neglected  since  the  opening  times  occur  muen 
faster  than  roll  damping  time  constants?  e.g.,  3.333  seconds 
vs.  3.33  seconds  respectively  for  a typical  situation. 
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Although  not  well  defined  in  magnitude,  it  is  possible 
that  the  retractei  fins  may  have  an  initial  coniltion  of 
opening  rate  due  to  an  impulse  transfer  either  from  the 
release  of  the  fin  latch  block  or  from  the  gun  barrel  gas 
dynamics.  Both  of  these  Influences  are  unknown  quantities 
that  may  occur  in  varying  degrees.  Another  concern  relates 
to  the  actual  fin  structural  dynamics,  and  a question  to  b? 
addressed  concerns  wnetaer  a fin  torsion  or  bending  mode 
might  be  significantly  excited  during  the  opening  and  hence 
provide  a magnification  of  fin  stresses.  Since  the  actual 
environment  of  the  E 3 IP  during  launch  may  be  conservatively 
described  as  “severe",  the  principal  source  of  flight  data 
to  answer  the  aforementioned  concerns  has  been  photographic 
in  nature. 

It  is  the  purpose  of  the  studies  described  herein  to 
develop  equations  of  motion  for  fin-opening  dynamics  under 
the  assumption  of  angular  momentum  conservation  such  that 
time  histories  may  be  obtained  for  the  opening  process  vith 
various  initial  conditions  of  roll  rate  and  fin  opening  rats 
for  both  rigid  and  elastic  fins.  From  computer  solutions  of 
these  equations  of  motion,  it  vould  be  hoped  that  some  con- 
clusions might  be  reached  relative  to  fin  stress  tims 
histories  in  order  to  provide  a tetter  insight  into  the 
overall  problem  that  is  loosely  categorized  under  the  title 
of  "fin-opening  dynamics". 
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II.  THEORY 


The  development  of  the  relations  describing  fin-opening 
dynamics  is  quite  easy  to  follow  when  one  employs  the  energy 
approaches  that  are  implicit  in  the  Lagrange's  equation 
technique;  e.g.,  Oolistein,  Ref.  1.  The  generalized  coord- 
inates vere  chosen  such  that  the  rigid  system  dynamics  were 
described  by  the  missile  rotational  angle,  <p  , ana  the  fin 
opening  angle,  © . Fin  elasticity  was  modeled  by  the  funda- 
mental bending  and  torsion  modes  using  a simple  pendulum 
analogy  where  bending  and  torsion  deformations  were  expres- 
sed by  the  coordinates  CL  and  0 respectively.  The  pendulum 
model  was  referenced  to  the  fin  base.  Figures  1 and  2 show 
the  coordinate  systems. 

Fin  geometry  is  illustrated  on  Fig.  3 in  accord  with 
Ref.  2.  Estimated  values  locating  the  center  of  gravity  of 
the  complete  fin  and  respective  mass  moments  of  inertia  were 
used  for  the  rigid  fin  model.  When  fin  elasticity  was 
modeled  (a  selectable  option),  the  mass  anl  inertial  proper- 
ties were  broken  down  to  considerations  of  fin  base  and  fin 
blade  elements  respectively.  Table  1 tabulates  the  values 
used  for  both  the  rigid  and  elastic  fin  analyses. 

One  may  conceptually  visualize  a refined  or  improved 
finite  element  elastic  model  using  the  fin  segmented  into  a 
four  degree  of  freedom  (d.o.f.)  system  as  sketched  in  Fig. 
4.  Aiding  more  degrees  of  freedom  will  improve  the  elastic 
modeling,  but  in  exchange,  the  equations  will  become  more 
complex.  Since  the  stiffness  values  for  the  pendulum  model 
arise  from  previous  dynamic  calculations  made  in  a station- 
ary frame  of  reference,  one  must  recognize  that  the  bending 
and  torsional  excitations  in  time  history  calculations 
represent  an  approximation  much  in  the  sense  of  using 
Talerkin's  method  in  aeroelastic  calculations  where  wini-off 
modes  are  used  to  develop  the  wind-on  modes  when  seeking 
in-flight  stability  boundaries. 


A.  Equations  of  Motion: 

In  terms  of  the  generalized  coordinates  4>  , © , oc  and 
P,  the  kinetic  and  potential  energies  of  the  missile  may  be 
established.  The  kinetic  and  potential  energy,  T and  V 
respectively,  may  be  functionally  recognized  as: 

r - r(  4>  , ©,oc , (3  , cj> , e ,6c,  (5,  t ) ...(  i ) 

and 

v « v(  4>  , e,(Y,  (3  , t)  . ..  ( 2 ) 
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TABLE  I 


Mass  and  Inertia  Properties 

ERSP  (.r/o  fins): 

I = 3.659  lb-in-sec2 

(estimated  about  missile  axis  of  symmetry) 

Fin  (one  of  six ) : 

a.  Mass:  MB  = 0.23  1437  lb-seca-inl 

MF  = 0.336134  lb-sec2 -in"1 
M = 0.33753:  lb-secz-in'  (total  value) 

b.  Location  of  c.g.  in  Fin  Coords.  (_2  = 2): 

Base : x 1 B = -3 . 142  in. 

i3B  = +3.301  In. 

Fin:  x ,F  = -0.316  in. 

1 3F  - +6.231  in. 

Total:  X)  = -0.774  in. 

x3  = +5.225  in. 

R0  = 2.330  in.  (fin  pi vot-to-ER9P  centerline) 

c.  Mass  Moments  of  Inertia: 

Total  fin  relative  to  total  fin  c.g.  location 
I , = 0.23742  lb-in-sec2 

12  = 0.39353  lb-in-sec2 

13  = 0.32533  lb-in-sec2 

Fin  (excluding  base)  relative  to  fin  c.g.  location 
I ,F  = 3.35329  lb-in-sec2 
I-gp  = 0.25321  lb-in-sec2 
I3f=  0.32492  lb-in-sec2 
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Fig.  I : Missile  Coordinate  System 


Fig. 2:  Elastic  Fin  Coordinates 
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Fig.  3:  Fin  Geometry 
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A:  two  BEND1NS  B:  two  torsion 

DEGREES  - Or-  FREEDOM  DEGREES  -OR-  FREEDOM 


Fig. 4 : Fin  Elasticity  Modeling 


For  the  freely  rotating  fissile,  the  roll  angle,  <£  , kill 
not  in  general  involve  potential  energy  storage.  The  fin 
opening  angle,  © , will  involve  potential  energy  storage 
upon  impact  Kith  the  crushabls  limit  stops.  The  fin  bending 
and  torsion  elastic  model  using  the  simple  pendulum  as  an 
approximation  to  the  fundamental  modes  will  introduce 
potential  energy  due  to  elastic  deformations  in  ol  and  (3  . 

The  application  of  Lagrange's  method  to  the  abov°  rela- 
tions for  kinetic  and  potential  energy  allows  us  to  express 
the  equations  as: 

_ iT 

for  q,  = <$>  , q2  = © 
q3  *c*  , q4  = $ 

The  above  equations  correspond  to  four  angular  momentum 
conservation  relations  with  coordinates  tailored  for  the 
problem.  5ince  some  of  the  coordinates  are  in  a moving 
frame  of  reference,  the  equations  include  terms  described  as 
Coriolis  type  accelerations.  Their  presence  is  automatical- 
ly included  in  the  Lagrangian  development.  In  the  theory  of 
small  vibrations,  the  derivative  of  kinetic  energy  with 
respect  to  the  generalized  coordinate,  £T/3qL,  is  omitted 
as  a consequence  of  linearization,  such  is  not  the  case  in 
the  nonlinear  relations  considered  here  since  the  q,  = <t>  and 
q2=  © deformations  were  not  small  perturbations.  The  right 
hand  side  of  eqn.  3 being  equal  to  zero  implies  that  the 
system  was  not  consiierei  with  an  external  forcing  function. 
However,  solutions  for  the  set  of  second  order,  nonlinear 
homogeneous  differential  equations  will  depend  upon  the 
prescribed  initial  conditions  i a q L a n i q , . 


Hl 


= 2 (i=l,  2,  3 and  4) 

...(  3 ) 


B . Solution  Procedure: 

The  terms  involved  in  eqn.  3 are  developed  in  Appendix 
A.  The  solution  procedure  for  numerical  analysis  involves 
expressing  the  equations  in  matrix  form  as: 

A q = g . . . ( 4 ) 

vhere  the  square  matrix  A introduces  the  inertial  properties 
of  the  system  with  a time  dependent,  nonlinear  coupling 
relative  to  the  fin-openiog  angle,  © , using  trigonometric 
relations.  The  matrix  A is  not  diagonal,  hence  convenient 
solutions  are  rot  readily  obvious.  The  column  matrix  g 
1 volves  tne  elastic  moments  due  to  fin  bending  and  torsion. 
It  also  involves  cross  product  terms  of  angular  rates  in  an 
interactive  manner  in  addition  to  a trigonometric  dependence 
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upon  fin  opening  angle. 

Nunerical  solutions  of  eqn.  4 were  obtained  using  tne 
IBM  363-6?  digital  computer  located  at  the  W.R.  Church  Com- 
putation Center,  Naval  Postgraduate  School.  The  solution 
procedure  was  based  upon  a fourth  order  Runge-Kutta  integra- 
tion procedure  using  the  features  of  the  Continuous  System 
Modeling  Program  (CSMP)  as  described  by  Speckhart  and  Sreen, 
Ref.  3.  Software  details  are  given  in  Appendices  B and  C. 

The  A matrix  in  eqn.  4,  although  time  dependent,  is 
nonsingular  at  any  given  instant  of  time  based  upon  physical 
arguments.  Therefore,  each  application  of  the  Runge-Kutta 
integration  scheme  required  uncoupling  the  highest  order 
time  derivatives  in  eqn.  4 by  performing  an  interim  step; 
i.e.,  at  any  time  t,  determine: 

..  -i 

q(t)  = A ( q , t ) g(q,  q,  t)  . . .(  5 ) 

Successive  applications  of  the  algorithms  available  by  the 
C5MP  procedure^proviied  straightforward  time  history  list- 
ings of  q and  q . An  output  sample  is  given  in  Appendlr  D. 

Since  the  CSMP  procedure  allows  the  use  of  FORTRAN  type 
logic  expressions,  it  was  relatively  simple  to  provide  a 
physical  modeling  of  the  crushing  action  for  a representa- 
tion of  the  fin  impacting  into  the  limit  travel  stops. 
Although  the  use  of  another  computational  language  may 
appear  foreboding  as  a usable  tool.  It  was  the  feeling  cf 
the  author  that  the  documentation.  Ref.  3,  and  principles 
involved  in  CSMP  were  so  well  expressed  th3t  anyone  moder- 
ately familiar  with  computational  methods  would  have  little 
difficulty  applying  the  technique. 


C . Fin  Holes  : 

The  desire  to  include  the  first  few  vibration  modes  of 
the  fin  into  the  time  history  dynamic  moiel  led  to  a search 
for  a finite  element  representation  of  the  fin.  Fortunat- 
ely, a finite  element  model  consisting  of  77  plate  elements 
with  92  node  points  was  available  to  represent  the  fin  as  a 
thin  plate,  and  the  required  information  was  provided  bi- 
Payne,  Ref.  4.  A planviev  s«etch  of  the  finite  element 
model  may  be  seen  on  Fig.  b.  The  inertia]  and  stiffness 
properties  of  the  fin  were  based  upon: 

p = Mass  density  - 7.323x13-4  1 b-see2 -in ^ 

E = Foung's  modulus  = 32,77xl2-<6  psi 

3 = Shear  modulus  = 11. 22x12+6  ps 1 
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Fig. 5:  Plate  Element  Layout  for  Fin 


The  fin  element  model  was  used  as  data  input  to  a structural 
analysis  program,  SAP  IV,  available  at  the  NPS  computer 
center  using  the  instructions  as  documented  by  Bathe,  Wilson 
and  Peterson,  Ref.  5.  The  first  three  modes  were  identified 
as  first  and  second  bending  and  first  torsion  at  character- 
istic frequencies  of  89.4,  411.4  and  517.0  Hz  respectively. 
Sketches  of  the  mode  shapes  may  be  seen  on  Figs.  6 to  3. 

The  fin  was  excited  on  a shaker  table  available  in  the 
Department  of  Aeronautics,  NPS  for  the  purpose  of  verifying 
the  modal  frequencies  and  modes.  The  actual  mounting  of  the 
fin  onto  the  table  was  by  means  of  clamping  the  fin  base 
between  two  - 1/2  inch  thick  aluminum  plates  such  that  the 
fin  was  harmonically  eicited  by  lateral  motion  of  the  base; 
i.e.,  oscillations  normal  to  the  fin  chord  plane.  •'Measured 
modal  frequencies  were  31.3,  355  and  430  Hz  with  the  modes 
identified  by  observing  the  nodal  lines  using  salt  crystals. 
The  ordering  of  the  mode  shapes  and  tie  location  of  the  node 
lines  were  in  accord  with  the  computer  results.  The  lo.°rPd 
values  of  the  characteristic  frequencies  as  compared  to  the 
theoretical  values  may  be  attributed  to  the  clamping 
mechanism  being  elastic. 

The  alternate  wedge  shaped  fin  (ftef.  3),  Fig.  9,  was 
also  vibrated  on  the  shaker  table.  Natural  frequencies  were 
experimentally  observed  at  84,  354  ani  422  Hz  for  first  and 
second  bending  and  first  torsion  modes  respectively.  The 
proximity  of  the  experimental  modal  frequencies  to  those  of 
the  diamond  shaped  fin  would  suggest  using  approximately  the 
same  theoretical  values  for  the  fin  dynamic  analysis. 

Equivalent  pendulum  model  stiffness  values  were  esti- 
mated in  the  fin-dynamic  analysis  based  upon  the  fin  inertia 
values  of  Table  I with  an  assumed  pendulum  pivot  point  in 
the  neighborhood  of  the  fin  base.  Based  upon  this  model, 
equivalent  spring  rates  for  approximating  the  first  bending 
ani  torsion  modes  were  (respectively): 

K = 1.581x13+4  in-lb/rai 

oc 

= 5.190x10+4  in-lb/rad 


The  first  and  second  bending  modes  could  be  approx- 
imated by  a double  pendulum  analogy  as  a further  refinement 
to  the  initial  analysis.  Based  up a a a knowledge  of  the  node 
line  for  the  second  bending  mode  ani  the  corresponding 
frequencies  for  first  ani  second  bending,  an  Intuitive  form 
of  parameter  identification  could  be  used  to  develop  the 
double  pendulum  model  with  respect  to  both  inertia  and 
spring  stiffness  traits.  This  type  of  a refinement  was  not 
attempted  in  the  dynamic  model  used  in  tne  following  studies 
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with  the  primary  reason  being  that  the  analysis  was  an 
exploratory  effort  to  evaluate  the  importance  of  fin  elasti- 
city upon  fin-opening  behavior.  The  initial  considerations 
of  only  the  fundamental  moles  should  be  viewed  as  a start- 
ing point  for  subsequent  analysis  refinements. 
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Fig.  6 . Fin  First  Bending  Mode 


REF.  CH'D.  PLfcNE 


[SECOWO  BESiD'M&] 

Moot  2’ 

All  .-4  WE 
uo^3  £56d  ^>e-C-1 


Fig. 7:  Fin  Second  Bending  Mode 


y = ia.ooo 
y * k.  a 20 

i ■ 0,0 


PIVOT  AvYVS 
COOCO.  SYSTEM 


V 


C.a.  LOCATION 


y ,r 


Fig.  8:  Fin  First  Torsion  Mode 
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Fig. 9:  Wedge  Fin  Geometry 
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III.  DISCUSSION 


A.  Background : 

The  fin-opening  dynamics  for  the  gun-launched  guided 
projectile  (ERSP)  occurs  after  the  spirning  EROP  leaves  the 
gun  barrel.  Upon  departing  the  barrel,  a slight  initial 
deceleration  of  the  vehicle  causes  the  release  of  a fin 
latch  or  lock  unit,  which  in  turn  releases  the  fins  ani 
permits  them  to  open  to  the  fully  deployed  position  under 
the  impetus  of  the  centrifugal  accelerations.  During  the 
deployment  process,  angular  momentum  may  be  considered  as 
being  conserved  since  other  dissipative  actions  such  as 
aerodynamic  roll  damping  are  not  of  major  significance 
during  the  initial  fin-opening  time  interval.  The  fins  open 
until  reaching  a 63  degree  position,  then  they  encounter 
small  cylindrical  stops.  The  fins  crush  the  stops  by 
plastic  deformation  ani  then  bio*  bacK  into  a detent  pos- 
ition corresponding  to  a 53  degree  position,  or  in  an  aero- 
dynamic sense,  a 30  degree  sweepback  angle. 

Details  of  the  fin-opening  process  nave  been  obscured 
because  of  inherent  measurement  difficulties  in  the  adverse 
experimental  environment  associated  with  launching  from  a 
Naval  rifle.  Photographic  observation  of  t h ? launch  process 
combined  with  inspections  of  the  ERSP  after  recovery  has 
provided  clues,  which  will  be  used  in  the  ensuing  analysis. 
One  significant  observation  from  fligat  tests  has  be  n that 
the  initial  roll  rats  of  the  ER3P  upon  departing  the  rifle 
was  dependent  upon  the  configuration  used  in  the  obturator 
rings.  The  obturator  rings  become  engrave!  into  the  rifling 
of  the  gun  barrel  during  launch.  The  relative  slippage 
between  the  rings  ani  the  ERJP,  which  is  controllable  to  a 
limited  extent,  acts  to  protect  the  ESJP  from  experiencing 
excessive  initial  roll  rates  at  launch.  However,  initial 
roll  rates  have  been  observed  up  to  33  revolutions  per 
second.  Basel  upon  these  observations,  initial  missile  roll 
rates  up  to  33  RP5  will  be  considered  during  the  parametric 
variations  of  the  analyses. 


B.  Crushable  Fin  Stop: 

The  rigid  fin  analysis  was  used  for  estimating  the 
magnitude  of  the  moments  induced  about  tae  fin  pivot  due  .o 
the  crushing  action  of  the  stops  during  fin  deployment.  A 
sketch  of  the  fin  stop  is  shown  on  Fig,  13.  In  estimating 
the  magnitude  of  the  moment  required  toticause  the  fin  stop 
to  become  crushed  into  a "washer-shaped"  object,  an  assump- 
tion of  a 123,333  psi  ultimate  compressive  stress  led  to  an 
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initial  guess  of  57,333  in-lb  total  moment  reaction  at  the 
fin  pivot  ails  for  the  sit  fins.  This  initial  value  of  the 
crushing-moment  level  provided  by  the  fin  stop  was  used  in  a 
parametric  study. 

Figure  11  summarizes  the  effects  of  tne  total  moment 
magnitude  (from  six  fins)  provided  by  the  crushable  stops  on 
the  maximum  value  of  opening  angle.  The  assumption  was  made 
that  fin  opening  began  with  an  initial  missile  roll  rate  of 
13  revolutions  per  second  and  zero  initial  fin  opening  angle 
rate.  Furthermore,  on  the  assumption  that  initial  contact 
with  the  fin  stop  occurred  at  a fin  opening  angle  of  63 
degrees,  the  fia  could  be  brought  to  a rest  condition  by  an 
additional  6 degrees  of  opening  angle  if  the  crushing  moment 
were  43,030  in-lb.  Several  things  to  note  in  the  analysis 
include : 

o The  original  guess  of  67,303  in-lbs  stopping  moment  was 
a startup  value  for  a parametric  stuiy  only. 

o Changes  in  initial  condition  of  missile  roll  rate  and 

fin  opening  angle  rate  would  require  different  values  of 
crushing  moment  if  the  constraint  of  bringing  the  flQ  to 
a rest  condition  in  5 to  13  degrees  more  opening  angle 
after  the  fin  contacted  the  stop  were  a valid  all 
encompassing  boundary  condition. 

o The  purpose  of  the  crushable  stop  logic  described  in 

section  II  was  solely  to  provide  a mechanism  for  bring- 
ing the  fin  opening  process  to  a rest  condition  in  a 
reasonable  manner. 

o The  elastic  fin  dynamics,  in  particular  the  fin's  elastic 
response  in  bending  and  torsion,  was  not  critical  during 
the  crushing  operation. 

Since  incorporating  the  logic  of  the  crushable  step  for 
the  purpose  of  arresting  the  fin  opening  angle  rate  was 
accomplished  in  a representative  manner,  and  did  not  lead  to 
any  critical  bending  or  torsion  responses,  further  investi- 
gation of  the  fin  moment  input  from  crushing  the  stop  did 
not  appear  warranted.  The  actual  initial  condition  range 
considered  in  the  investigations  varied  in  the  following 
manner : 

o t = 3 missile  roll  rate  varied  from  3 to  33  RPS 

o 6(3) ; t=3  fin  opening  rate  varied  from  2 to  43  rad/sec 
(2292  ieg/sec ) . 

The  effects  of  the  above  range  of  initial  conditions  upon 
the  crushing  of  the  fin  stop  and  consequently  the  maximum 
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opening  angles  are  shown  on  Fig.  12  for  tne  situation  of 
rigii  fin  dynamics.  As  will  be  noted,  considering  the  stop- 
ping  moment  from  crushing  action  as  constant  (42x13+3  in-lb) 
did  not  prevent  the  maxiumum  opening  angle  from  exceeding  92 
degrees  when  the  initial  roll  re'.e  exceeded  32  RPS.  A more 
refined  analysis  might  have  modeled  the  stopping  moment  due 
to  plastic  crushing  of  the  stop  as  being  a gradual  buildup 
in  moment  value.  The  computer  logic  to  implement  such  a 
feature  is  not  difficult,  but  was  not  employed  since  the 
maximum  fin  elastic  deformations  in  general  appeared  before 
the  fin  stop  was  contacted  during  deployment. 
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Fig.  10’  Sketch  of  Fin  Stop 
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C . Opening  Angle  Rate : 

The  time  in  milliseconds  for  the  fins  to  reach  an 
opening  angle  of  53  degrees  is  shown  on  Fig.  13  as  a 
function  of  initial  opening  rate  for  various  initial  values 
of  missile  roll  rate.  Experience  from  test  firings  of  the 
ERJP  has  shown  that  the  condition  of  the  interface  between 
the  obturator  rings  and  the  missile  body  allows  a reasonable 
degree  of  control  upon  the  initial  missile  roll  rate.  The 
opening  angle  rate,  which  can  be  noted  on  Fig.  13  as  having 
a strong  influence  upon  the  fin-opening  times,  is  not  as 
readily  amenable  to  control  during  testing.  Possible 
factors  influencing  the  intial  fin-opening  rate  would 
certainly  include  muzzle  gas  dynamic  interference  effects  as 
well  as  impulsive  loadings  induced  by  release  of  the  fin 
latch. 

The  latter  effect  has  been  estimated  in  termst>  of  the 
impulse  created  by  a latch  delta-velocity  (AV)  "jump”  o' 
one  foot  per  second  during  latch  release.  The  analysis, 
details  of  which  may  be  found  In  Appendix  E,  would  Indicate 
a sensitivity  of  0.66  rad/sec  from  an  one  foot/sec  velocity 
"jump”  at  release  for  a latch  weighing  approximately  3.51 
pounds. 

Unfortunately,  it  is  difficult  to  estimate  the  transfer 
of  impulse  from  the  latch  into  fin-opening  rate.  However, 
there  is  experimental  evidence  supporting  the  notion  that 
the  fin  latch  moves  forward  relative  to  the  ER3P  upon  the 
system  leaving  the  gun  barrel.  The  determination  of  what 
fraction  of  the  latch's  relative  forward  motion  is  used  as 
an  impulse  transfer  to  initialize  the  fin-opening  rate  is  a 
difficult  item  to  quantify  since  the  test  environment  in  the 
presence  of  the  gun  muzzle  blast  is  quite  advTse.  From  an 
order  of  magnitude  viewpoint,  a velocity  "jump"  from  the 
latch  of  53  feet  per  second  for  impulse  transfer  would  seem 
appropriate  considering  the  fact  that  the  latch  can  er, count- 
er an  accelerating  force  due  to  the  gun  barrel  gas  flow 
(relative  to  the  missile)  upon  launch. 

The  inertial  proper tlesitof  the  fin  were  revised  to 
reflect  the  alternate  "wedge”  fin  shape,  Ref.  6,  and  a rigid 
fir  opening  dynamics  analysis  was  performed.  The  effect  of 
using  the  wedge  fins  upon  the  time  to  reach  a 63  degree 
opening  angle  is  shown  on  Fig.  13  by  the  dashed  curves.  The 
wedge  fins  generally  required  slightly  more  time  to  reach  63 
degrees  relative  tc  the  diamond  cross-section  fin.  Ref.  2 , 
but  in  general  the  time  difference  vas  on  the  order  of  a 
five  percent  or  less  increase. 
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D.  Elastic  System  Response: 


A typical  time  history  plot  showing  fin  opening  angle 
for  various  opening  rate  initial  conditions  at  an  initial 
missile  roll  rate  of  13  revolutions  per  second  (RP’S)  is 
shown  on  Fig.  14  for  a typical  situation  of  an  elastic  fin 
response.  The  contact  by  the  fin  with  tie  crushable  stop  at 
63  degrees  opening  angle  is  followed  by  a short  period  of 
deceleration  to  a zero  opening  rate  condition  and  is  denoted 
on  Fig.  14  by  the  dashed  fairing.  Fin  elastic  effects  were 
not  noticable  upon  the  curve  shape  since  angular  momentum  in 
the  theta  <©)  direction  was  not  directly  influenced  by  fin 
bending  and  torsion  response.  The  time  to  reach  60  degrees 
of  fin  opening  was  increased  slightly  when  fin  elasticity 
was  included  in  the  dynamic  modeling  as  a result  of  fin 
bending  effects  naving  an  influence  upon  missile  roll  rate 
time  histories,  which  in  turn  interacted  with  fin  deployment 
through  a centrifugal  force  dependence.  A first  order  con- 
sideration of  opening  angle  time  histories  may  be  viewed  as 
being  reasonably  well  represented  by  the  rigid  system  model 
since  fin  elasticity  is  not  a strong  influence  upon  the 
theta  (©)  generalized  coordinate. 


A significant  factor  in  fin  deployment  is  the  mass 
moment  of  inertia  increase  about  the  roll  axis.  Deploying 
the  fins  to  the  53  degree  position  typically  increased  the 
9-lnch  ERGP  roll  axis  moment  of  inertia  by  25  percent. 
Angular  momentum  considerations  would  lead  one  to  expect  a 
22  percent  decrease  (3.53  = 1/1.25)  In  missile  roll  rate 
relative  to  the  initial  value  at  departure  from  the  gun 
barrel. 


The  rigid  fin  roll-rate  time  history  plot  is  shown  on 
Fig.  15  for  an  initial  roll  rate  of  13  RPS  and  varying  con- 
ditions of  initial  opening-angle  rate.  Although  the  times 
to  reach  the  60  degree  deployment  situation  were  dependent 
upon  initial  conditions,  the  ratio  of  deployed  system  to 
initial  system  roll  rates  remained  invariant  by  the  angular 
momentum  conservation  considerations  stated  above.  It 
should  be  noted  that  these  remarks  are  applicable  only  to 
the  system  with  the  rigid  fin  assumption. 

Figure  16  covers  the  same  range  of  initial  conditions 
as  in  the  preceding  figure,  except  that  the  dynamic  terms 
from  the  fin  motion  included  an  elastic  modeling  influence 
from  first  bending  and  torsion  modes  respectively.  The  time 
history  curves  have  the  time-to-reach  53  degrees  flagged 
out,  and  comparison  between  Figs.  15  and  15  will  substan- 
tiate the  earlier  comment  that  fin  elasticity  effects  were 
not  a dominant  influence  upon  fin  opening.  But  Fig-  15  does 
show  that  the  inclusion  of  fin  elasticity  significantly 
altered  the  roll  rate  time  history,  and  in  particular,  the 


value  of  t) i s s 1 1 s roll  rate  when  the  fins  were  deployed  to  50 
degrees  was  not  invariant  with  fin-opening  rate  as  was  the 
case  for  the  rigid-fin  dynamic  model.  The  influence  of  fin 
elasticity  may  be  attributed  to  a coupling  between  the  fin's 
bending  mode  and  the  missile's  roll  mode.  Fin  bending 
motion  can  be  visualized  as  interacting  with  angular  momen- 
tum about  the  roll  ails  with  an  added  dependence  from  the 
fin  deployment  angle. 

These  remarks,  concerning  the  influence  of  fin  elas- 
ticity upon  roll-rate  time  histories  suggest  that  diffi- 
culties might  be  experienced  in  estimating  fin-opening  per- 
formance by  photographic  observations.  It  would  be  quite 
possible  to  mistake  the  influence  of  fin  bending  excitation 
as  an  unknown  error  (or  "noise")  in  missile  roll  rate.  In 
addition,  actual  photographic  evidence  of  missile  roll  angle 
(and  hence  rate)  may  occasionally  be  clouded  by  the  gun 
barrel  blast.  It  might  be  difficult  to  correlate  an  obtur- 
ator ring  configuration  with  a transient  roll-rate  value, 
expecially  if  the  missile  roll  rate  were  to  change  by  20  to 
30  percent  during  the  typical  23  to  30  millisecond  opening 
process . 

Another  interesting  consideration  is  provided  by  Fig. 
17,  which  shows  the  dependence  of  missile  Initial  condition 
upon  roll  angle  when  50  degree  fin  deployment  angle  is  first 
attained.  The  family  of  roll  angle  curves  would  Indicate 
that  as  a general  rule,  the  fins  reach  the  deployed  position 
in  less  than  one-third  of  a revolution.  One  point  to  note 
is  that  the  roll  angle  for  50  degrees  fin  deployment  was 
invariant  with  missile  initial  roll  rate  providing  that  the 
fin-opening  rate  was  initially  zero.  However,  this  latter 
observation  may  be  somewhat  academic  since  the  evidence  from 
test  observations  during  firing  tends  to  support  an  assump- 
tion of  a non-zero  initial  fin-opening  rate. 

Representative  elastic  fin  bending  response  time 
histories  for  an  initial  missile  roll  rate  of  10  RP5  may  be 
seen  on  Fig.  18  for  various  opening  rate  initial  conditions. 
Of  particular  note  is  that  the  bending  response  apppars  as  a 
superposition  of  a time  varying  waveform  upon  a linearly 
Increasing  variation  with  the  amplitude  of  the  wave  being 
dependent  upon  initial  opening  angle  rate.  As  may  be  noted, 
peak  values  of  bending  deflection  occur  on  the  second  wave 
for  the  lower  values  of  initial  opening  rate  and  on  the 
first  wave  at  the  higher  values  of  Initial  condition.  At 
the  larger  values  of  initial  opening  rate,  the  time  to  reacn 
50  degrees  of  opening  angle  was  20  milliseconds  or  less  and 
hence  the  response  could  accommodate  only  one  wave  peak 
before  the  action  of  the  crushable  fin  stop  entered  into  the 
dynamics . 
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Although  the  perioi  of  the  wave  was  not  clear  on  the 
curves,  Fig.  18,  an  approximate  estimate  of  1?  milliseconds 
for  a wave  perioi  would  indicate  the  presence  of  a 59  Hz 
frequency  in  the  transient  response.  The  fundamental 
frequency  of  the  fin  in  a non-rotating  coordinate  system  was 
estimated  as  39. d Hz  (c.f.,  section  II-C),  which  serves  to 
illustrate  that  the  rotating  coordinate  system  does  have  an 
influence  upon  the  modal  frequencies  and  consequently,  also 
upon  the  transient  responses.  Further  investigations  of  the 
transient  response  behavior  certainly  should  consiier  the 
effects  of  higher  bending  modes  as  well  as  a sensitivity 
study  due  to  changes  in  the  various  stiffness  constants. 

The  transient  torsional  response  for  the  sane  range  of 
initial  conditions  as  described  above  may  be  found  on  Fig. 
19.  The  curve  shapes  are  similar  to  those  shown  for  bending 
response!  however,  the  response  magnitude  appears  reduced  by 
three  orders  of  magnitude.  This  observation  would  lead  one 
to  consider  dynamic  modeling  of  fin  elasticity  as  being 
primarily  dominated  by  fin  bending  response. 

The  main  impetus  tor  including  fin  elasticity  in  the 
dynamic  modal  of  fin  opening  was  to  identify  situations  when 
fin  root  bending  stresses  could  become  sufficiently  large  so 
as  to  induce  structural  yielding.  Figure  22  summarizes  the 
peak  values  of  the  generalized  first  bending  coordinate  for 
the  complete  range  of  initial  conditions.  An  alternate 
scale  is  also  presented  of  maximum  induced  tending  stress 
encountered  at  the  fin  root.  The  corresponding  maximum 
stress  scale  was  based  upon  translating  benilng  deflection 
to  bending  moment  using  the  effective  spring  stiffness 
constant,  , described  in  section  II-C.  If  yielding  were 
assumed  to  occur  in  the  neighborhood  of  152,000  psl  tensile 
stress,  tnen  evidence  of  permanent  set  in  the  fis  bending 
coordinate  could  be  expected  if  the  Initial  missile  roll 
rate  ezceedei  15  to  22  RPS.  It  should  also  be  noted  that 
the  maximum  stress  encouuter  was  influenced  more  strongly  by 
initial  roll  rate  rather  than  by  initial  fin-opening  rate. 

In  order  to  emphasize  the  latter  comment.  Fig.  23  has  been 
prepared  to  illustrate  the  peak  flu-bending  induced  stress 
encounter  as  a function  of  Initial  missile  roll  rate  for 
several  cases  of  initial  fin-opening  rate. 

Figures  22  and  21  illustrate  that  the  most  effective 
way  of  avoiding  missile  residual  rail  rate  due  to  permanent 
set  of  the  fins  Is  to  configure  the  obturator  rings  so  tba-. 
initial  missile  roll  rates  are  kept  below  approximately  15 
RPS. 
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©Co)  sec 

Fig.  20:  Peak  Fin  Bending  vs. 
Initial  Opening  Rate 


A. s an  aside,  it  should  be  recognized  that  the  linear 
stress-strain  relationship  of  the  fin  material,  which  was  an 
inherent  feature  of  the  fin  elasticity  model,  vouli  not  be  a 
valid  assumption  for  higher  values  of  initial  conditions 
where  material  yielding  night  be  expected.  However,  the 
curves  of  peak-stress  encounter  represent  events,  which  for 
stresses  beyond  yielding,  only  occur  for  a few  milliseconds . 
Therefore,  improvements  in  the  dynamic  modeling  should  also 
reflect  material  property  dynamics  lx  one  were  seeking  to 
quantify  estimates  of  fin  yielding.  Fortunately,  the  CS^P 
technique.  Ref.  3,  allows  the  convenient  inclusion  of  logic 
statements  that  can  accommodate  improved  dynamic  modeling  in 
case  the  scope  of  the  investigations  were  broadened. 
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E.  Experimental  Correlation: 


As  one  would  expect  from  earlier  comments  male  in  this 
report,  the  acquisition  of  experimental  information  for  des- 
cribing fin-opening  iynamics  was  potentially  limited,  both 
in  amount  and  scope,  because  of  the  adverse  nature  of  the 
launch  environment.  A limited  amount  of  fin -opening  perfor- 
mance has  become  available  based  upon  photographic  records. 
Results  of  this  technique,  described  in  Ref.  ?,  showing  the 
time  for  52  degree  fin  deployment  are  summarized  on  Fig.  22, 
It  should  be  noted  that  these  results  are  relative  to  the 
wedge-shaped  fins  which  are  slightly  different  in  properties 
from  the  diamond-shaped  fins.  The  latter  fin  configuration 
was  the  principal  consideration  in  this  report. 

The  two  faired  curves  on  Fig.  22  represent  results  from 
fin-opening  spin  tests  conducted  under  controlled  laboratory 
conditions  and  from  gun-launched  test  firings  using  an 
9-inch  rifle  available  at  the  N*?C  gun  range.  The  spin  test 
results  were  assumed  to  represent  a zgro  initial  condition 
situation  for  fin-opening  rate,  while  the  fin-opening  rate 
was  an  unknown  for  the  actual  gun-lauQchei  tests.  Table  II 
compares  predicted  results  from  the  dynamic  modeling  of  a 
(rigid)  wedge-shaped  fin  with  the  data  points  on  Fig.  22. 

For  the  situation  of  the  spin  tests,  the  predicted  values  of 
Tt0  were  in  general  slightly  lower  by  a few  milliseconds, 
using  the  assumption  that  the  initial  opening  rate  was  zero. 

The  actual  gun-launched  tests  provided  values  of  Tto 
that  could  be  used  in  an  indirect  manner  to  obtain  estimates 
of  in-sltu  initial  fin-opening  rates.  By  matching  the 
measured  and  predicted  fin -opening  times,  initial  fin- 
opening  rates  were  deduced  in  the  range  of  33  to  53  rads  per 
second;  c.f.,  Table  II.  These  results  add  support  to  the 
intuitive  feeling  that  the  opening  rate  was  due  to  an 
effective  fin  latch  av  "jump"  of  approximately  53  feetMper 
second,  where  it  is  recognized  that  the  term  "effective"  was 
used  to  cover  momentum  input  sources  to  fin  opening  from  gun 
blast  dynamic  Influences  as  well  as  the  kinematics  of  the 
fin  latch  mechanism.  In  any  event,  there  does  appear  to  be 
credence  provided  relative  to  the  range  of  fin-opening  rates 
considered  in  the  analysis  of  the  report,  as  veil  as  to  the 
results  of  the  analysis. 
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Fig.  22  Experimental  Fin-Opening  Times 
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IV.  CO'JCLUDINS  RE'URKS 


An  analysis  of  fin-opening  dynamics  has  been  completed 
for  the  9~incti  ERJP  using  four  generalized  coordinates  in 
conjunction  with  Lagrange's  equations  of  motion.  Phis 
method  procedurally  simplified  the  derivation  for  the 
coupled  equations  of  motion  into  an  algebraic  manipulation. 
The  elastic  model  for  representing  fin  dynamics  included 
both  fin  .bending  and  torsion  response  representations. 
Analyses  disclosed  that  the  bending  responses  dominated  the 
elasticity  influences  by  several  orders  of  magnitude,  which 
result  was  in  accord  with  Independent  analyses  conducted  by 
Boone,  Ref-  7,  cr  a 5-inch  gun-launched  projectile  with  a 
similar  fin  configuration. 

ihe  question  of  result  credibility  was  addressed  by 
making  a comparison  with  limited  experimental  results  giving 
opening  times,  and  from  these  comparisons  an  aided  benefit 
was  ottainedi  namely,  on  Indication  that  the  initial  fin- 
opening  rate  was  in  the  neighborhood  of  35  to  53  radians  per 
second  . 


Since  the  results  of  the  report  appear  physically  con- 
sistent, one  may  conclude  that  the  analysis  method  does  pro- 
vide tne  systems  designer  with  a tool  to  describe  fin- 
opening  dynamics,  including  Other  configurations,  in  a 
straightforward  manner.  dinisight  obtained  from  the  results 
would  suggest  that  improvements  Id  the  elastic  model  could 
be  made  by  adding  more  fin  bending  degrees  of  freedom  and  by 
incorporating  improved  logic  to  accommodate  the  effects  of 
material  yielding. 


The  residual  roll  rate  for  the  3-inch  ERdP  following 
the  completion  of  the  gun  launch,  fin-opening  transients  and 
roll  rate  decay  by  aerodynamic  damping  is  governed  by  both 
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structural  yielding  due  to  launch  and  opening  dynamics.  The 
latter  effects,  described  as  a main  goal  of  these  investiga- 
tions, appeared  to  be  dependent  primarily  upon  the  Initial 
missile  roll  rate  and  secondarily  upon  the  initial  condition 
of  fin-opening  argls  rate;  c.f..  Fig.  2'.  It  would  appear 
that  most  of  the  problems  associated  uitli  inertially  induced 
structural  yielding  by  the  fins  could  be  avoided  if  the 
obturator  ring  configuration  between  the  missile  and  the  gun 
barrel  were  defined  for  initial  missile  roll  rates  below  15 


RP3.  This  philosophy  does  not  eliminate  the  influence  of 
gun  barrel  gas  dynamics  and  fin  latch  momentum  transfer  into 
the  fin's  initial  opening  rate,  but  rather  it  relegates 
these  problems  into  a subordinate  role,  which  certainly  is 
to  be  desired  when  trying  to  achieve  engineering  design 
success . 


APPENDIX  A 


Equations  of  Motion: 

The  Lagrange  equation  approach  will  be  use!  for  setting 
up  tne  governing  relationships  for  the  fin-dynamics  problem. 
The  equations  will  include  terms  to  provide  an  elementary 
treatment  of  the  elastic  fin  response;  however,  it  will  be- 
come evident  in  the  analysis  as  to  which  terms  need  to  be 
removed  for  reversion  to  a rigia-fin  situation.  The  first 
step  in  the  procedure  is  to  locate  the  center-of-gravi ty 
(c.g.)  positions  of  the  various  fin  elements  relative  to  an 
inertial ly  oriented  reference  frame.  The  next  step  relates 
to  the  development  of  the  system  kinetic  energy  terms,  com- 
posed of  translational  ani  rotational  components  relative  to 
the  individual  c.g.  locations.  The  thiri  step  involves  the 
establishment  of  potential  energy  type  terms.  Finally,  the 
equations  of  motion  become  lientifiei  by  direct  substitution 
of  the  appropriate  derivative  terms  in  Lagrange's  equation. 


The  contiibution  to  energy  from  the  rotating  missile 
body  (without  fins)  is  given  by: 


AT  s 'a  T ,4,  2 

«— i * ’ ^ -L 

and 

A V = 0 


. . . (k.l) 


In  treating  the  fin  trauslatiooal  contribution  to 
energy,  assumptions  made  will  include: 


o Fin  is  modeled  by  a rigid  base  plus  an  elastic  fin  blade 
section.  The  fin  base  is  pin  connected  to  the  missile 
body  at  a distance  R„  from  the  missile  axis  of  symmetry. 

o Fin  elasticity  considers  only  the  fundamental  bending  ani 
torsion  modes  using  a simple  pendulum  model. 

o Fin  base  dynamics  can  be  obtained  from  rigid-fin  modeling 
by  appropriate  revisions  to  the  c.g.  ani  mass  terms. 


Fin  Base  Kinetic  Energy  

The  absolute  position  of  the  fin  bass  o.g.  relative  to 
the  inertial  coordinate  frame  of  reference,  x'  , and  x8 
(Fig.  1)  for  the  base  may  be  expressed  as: 


R0+ xiacos © + x3b S'N®  cos 
R0+  X,8COS©  + X3B=,\tvi©J-5lN4> 


...(A. 2) 


*3©  ' -X,aSlNG  + X3bCOS9 
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where  numerical  values  for  R0  , x,B  and  x3B  may  be  found  in 
Table  I.  The  next  step  is  to  take  the  time  derivative  of 
the  inertial  position,  followed  by  summing  the  squares  to 
obtain  the  square  of  the  velocity;  i.e., 

vf  = (x'3  ) + (x  '2b  f + (i  3B  )2  ...(A.  3) 

The  translational  contribution  to  kinetic  energy  for 
six  (5)  fin  bases,  following  algebraic  simplification, 
becomes:  . 

ATbt«c  6/2  Mb  VB^ 

C * 2r  2 2 2 ^ 2 

•'  = 3MB|c}>|Ro4-  X1b  COS  © 4-  X3b  sin  © + 

+ 2 Ro(x,b  COS  © -V  X36S1MG»')  4-  X|BX30SllM  2©] 

+-©2[x.l  + X32e]]  ...(A.4) 

The  rotational  contribution  to  kinetic  energy  requires 
expressing  the  angular  velocities  of  the  fin's  x,  , x2  and  x3 
coordinate  system  ( u_>, , <jJ2andu>3  respectively)  in  terms  of 
an  inertial  frame  of  reference.  Two  coordinate  rotations 
are  involved  to  describe  the  rigid  fin  angular  rotation  con- 
sisting of  first,  a rotation  about  tas  missile  axis  of 
symmetry  followed  by  a rotation  & about  the  fin  pivot  axis. 
It  should  be  recognized  that  this  is  considerably  simpler 
than  the  Euler  angle  approach  which  is  the  classical  treat- 
ment applied  in  the  more  general  situation  when  developing 
angular  momentum  conservation  relations  for  a rigid  body,  be 
it  either  a gyroscope  or  an  airplane. 

Consequently  one  obtains  that: 

« 

UJ,  = - (p  sin  o 

bO  2 — © . . . ( A . 5 ) 

u)3  = <$>  COS© 

which  provides  the  entries  for  the  rotational  contributions 
to  kinetic  energy  (for  5 fins)  as  follows: 

ATBro=  3[lie('Ct>  S\N©f+  I2B©2+I3e<<i>COS©2]  ...(A. 5) 

Fin  Base  Potential  Energy  

The  contribution  of  the  fin  base  model  to  system  poten- 
tial energy  would  depend  upon  the  relative  importance  of 
terms  such  as  the  eartn's  gravitational  field  or  fin  elastic 
deformations.  The  former  was  considered  as  a negligible 
factor  for  the  rotating  missile  environment  while  the  latter 
effect  was  precluded  by  tae  1 ac tt  of  identifiable  springs  in 
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the  base  area.  One  shouli  note  that  there  is  a orushabie 
stop  present  in  the  system  vhicli  absorbs  energy  after  the 
fin  opening  angle,  6>  , has  reached  62  degrees.  The  logic 
required  for  this  effect  will  be  treated  separately  later. 
Therefore,  based  upon  these  arguments,  we  shall  assume  that: 

A VB  =3  ... (A,?) 


Fin  Blade  Kinetic  Energy  

In  modeling  fin  elasticity  according  to  the  stated 
assumptions,  the  fin  was  considered  to  deflect  in  bending 
deformation  as  an  equivalent  pendulum  with  angle  oc.  about  a 
virtual  pivot  locatei  on  the  chord  plana  at  the  top  of  the 
fin  base.  Symmetry  conditions  were  assumed  for  the  elastic 
fin  about  the  miichord  line.  Hence,  the  absolute  position 
of  the  elastic  fin  c.g.  relative  to  the  inertial  coordinate 
frame  of  reference,  x,  , x2  and  X3  (Fig.  1)  may  be 
expressed  in  a form  similar  to  eqn.  A. 2 as: 


x 


IF 


X 


2F 


Ro+Xip  cos©  + X3f  S\M©JCOS$4AX3S\N  <J>  CX 
R0+XiP  COS©  4-  V3FSVN©JsiH(|)-AX3COS^Ob 


x'3f  = -X,rSlN0  + X3(r  COSO  ...(A. 3) 

The  numerical  values  involved  in  the  above  relations 
are  stated  in  Table  I.  Tne  term  Ax3  is  the  distance  in 
the  x 3 direction  (spanwise  for  the  fin)  from  the  virtual 
pivot  of  tne  simple  psniulum  model  to  tne  elastic  fin  c.g.. 
The  above  equations  include  a small  angle  assumption  upon 
the  bending  angle  deformation?  i.e.,  OL  = sinCX  . 


The  translation  contribution  of  the  six  (6)  elastic 
fins  to  kinetic  energy  may  be  evaluated  in  a manner  similar 
to  the  development  that  led  to  eqn.  A. 4.  Results  were: 

ATFtr=  3 {4>2[r^4  X,fCOS2©  4 x|F  S\N*©  4 

+ R0 (xlpCOS  © 4 X3f  SIM  &)  4XifX3fS'M  2©4Axfo£2J 
4 ©2[x,p  4 X.;]  + A X3  OC2 

+ 2<£©[-x,fs\m©  4 x3f cose] £1X3 cx) 

- 206c[eo4  XlpCOS  ©4  X3FS\Me](AX3)|  * * 

For  the  elastic  fin  model,  the  angular  velocity  com- 
ponents of  eqn.  A. 5 become  modified  by  inserting  bending  and 
torsion  angle  deflection  rates,  as  noted  below: 
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u),  = 


- 45  S\N  © 4 (k 
(jJ2=  © ...(A.  12) 

u)3=  <p  COS©  ■+  ^ 

The  contribution  of  rotation  to  Kinetic  energy  for  the 
six  elastic  fins  is  obtained  by  applying  eqn.  A. 12,  which 
yields : 

AT^  3 f Ilp  (-cj>  S inj  © + oc)2  + I2(r  ©2  + 

+ l3p(  4>  COS  © + (i)2  J ...(A. 11) 

Fin  Blade  Potential  Energy  — - 

The  elastic  modeling  of  the  fin  as  a sinple  pendulum  in 
both  bending  and  torsion  alius  us  to  express  the  potential 
energy  for  the  six  fins  as: 

AVF  = 3 [ C*2  4-  £2]  ...(A. 12) 

In  the  above  expressions,  the  estinate  for  the  elastic 
spring  constants  were  obtained  from  tne  separate  free  vibra 
tion  calculations  of  the  finite  element  fin  model  using  the 
structural  analysis  program  Known  as  SAP  IV. 
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The  last  iten  for  inclusion  ir  the  energy  fornulation  is 
the  crushable  stop  which  is  used  to  bring  the  fin  to  a rest 
condition  during  the  opening  process.  Tne  crushable  stop, 
one  per  blade,  is  a snail  steel  cylinder  positioned  on  tne 
nissile  base  frame  in  a manner  such  that  the  fin  base  makes 
initial  contact  at  a fin  opening  angle,  © , of  approx- 
imately 53  degrees  (1.217  radiaas).  The  fin  contact  with 
the  stop  is  initially  an  elastic  deformation  over  a small 
change  in  opening  angle,  but  once  material  yielding  has 
occurred,  tne  influence  would  correspond  to  plastic  deforma- 
tion. Therefore,  the  vorK  from  the  absorption  term  will  be 
modeled  as  an  equivalent  potential  energy  term  by  an  expres- 
sion for  the  six  fins  as  stated  below: 


AVS  = 


for 


© .LT.  1.24?  rad 


6 Ks(e-  1.347)  for  © .5E.  1.247  rad 
ani  iO/at  .GI.  2 


. ..(A. 13) 


where  Ks  is  an  estimate  of  the  stopping  moment  developed  by 
the  plastic  type  deformation  of  the  crushei  stop  during  the 
deceleration  of  the  fin  opening  angle.  The  actual  digital 
computation  will  cease  (by  a logic  test)  wnen  the  angular 
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rate  of  opening  angle  reaches  zero.  In  practice,  the  fins 
will  then  blow  back  due  to  the  aerodynamic  type  moments  into 
ttie  detent  (loosed)  position  of  60  degrees. 


Summary  of  Energy  Terms  

The  virtues  of  a systematic  energy  formulation  for  a 
complex  system  tv i 1 1 now  become  evident  since  the  establish- 
ment of  the  appropriate  coefficients  in  the  equations  of 
motion  becomes  simplified  to  a routine  algebraic  accounting 
procedure.  First,  the  kinetic  energy  terms  may  be  summar- 
ized in  their  entirety  as: 

T =AT  + A?Btr  - ATBrot  + ATFtr  « ATFrot  ...(A.U) 

wnere  the  contributions  to  the  right  nand  side  are  defiosd 
by  eqns.  A..1,  A. 4,  A. 6,  A. 9 and  A .11  respectively. 


Next , 
V 


trie  potential 


energy  terms  in 
+ AVr  + AV  e 


tneir 


entirety  are: 

. . . (A. 15) 


where  the  contributions  to  the  right  nand  side  have  been 
established  by  eqns.  A.l,  A.?,  A. 12  and  A. 13  respectively. 


Definition  of  Coefficients  

Equation  3,  whicn  is  Lagrange's  equation  for  establish- 
ing tne  differential  equations  that  describe  the  fin-opening 
dynamics,  is  restate!  below  for  sake  of  clarity: 


d_/ }J_\  r) T 

citV^  / aq. 

where 


(i 


1,  2,  3 and  4) 
. . . ( 3 ) 


q,  * 4>  . q?  = © 

q3  = a , q«  = 8 


The  above  set  of  four  coupled  differential  equations 
may  be  expanded  term  ty  term  using  the  kinetic  and  potential 
energy  expressions  of  eqns.  A. 14 

eqn . 


of 

of  the  first  energy  term  of 
a matrix  format  as  fellows: 

±(£L) 

dt 

where  the  square  matrix  A will 
con-diagonal  in  character.  Perm 
kinetic  energy  expression  yielded 
for  the  A matrix: 


and  A. 15.  In  the  expansion 
3,  it  is  convenient  to  use 


. ..(A. 16) 

in  general  be  symmetric  but 
by  term  expansions  of  the 
the  following  coefficients 
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A 2c ~ & ("  lg  +■  ^ ( X i*"  'll 

A 33-  6 [It(r  + 

A 4^-  £>  l3jr 

A (2  = A2i  * G MF{  X3f  COS  © -- X;k.-SIN  ©jA  X3  Ckf 

A ,3=  A3,  = - ^fMpR0AX3+[llr+NAFX3f  AX3]sm©  + MP X,FAX3COS©J' 

A \4  = A 4i  = 6>  I3f  COS  & 

A 23'  A32  “ 0 
A 24  = A 42  = <3 

^ ■*  A43'  ^ •»«(Avlo^) 

It  is  convenient  to  rearrange  tne  equations  of  notion  into 
matrix  forn  (eqn.  4)  as: 

N fi)  * {*)  ••■(4> 

where  the  column  vector,  g , is  derive!  from: 

Tq  “ ==^>  (b'l  = is(lK  q»  t>|  • . • ( a . 1 ? ) 

The  first  component  of  tne  g column  vector,  when  expanded, 
becomes : 

S'  * >»wi,V  bra^,q5+  t'jj'jj  + b?3g233  + bMq2g4 
where  ...(A. 17a) 

b ,2  = - & {[(T,-I3)  + M(  x|-  xf  )]siN  2©  4?M  R0[x3cos  © - X,  sm  ©]  + 

+ ?[M'Bxi8VMfV^f]co5  2«j 

bi 3 * -12  Mp  Axf  a 

G NApj^JvjpCOS©  4-  X3f  SIN  ©J  A X3  £X 
b23=  G IVp  COS  © 

°24-  &I3fSIN0 

The  second  component  of  the  g column  vector  is: 

c^t  + c^<i(<53+  9^  + -f*.  ...(A. 17b) 
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c(i  = 3 [(It- I3H  M(x|- X,2)]sin  2©  + MR0[x3cose-X,SiN@]  + 

4-  <®[mbXisX3b*'-MfXifX3f]COS 
cl3  = -6>IlFCOS© 

CIA  = ~<o  I3fs\m  © 

2 for  q,,  . I,T.  1.247  rad 

f*  - ■ 

-6  Ks  for  q2  .JE.  1.247  rad  and  dq2/dt  ,5T.  2 
rue  third  component  of  the  g column  factor  is: 

£3  = 4^  ’+<^12  ^ "*■  "^33  ...(A. 1?C) 

where 

d„  * <o  MPAX32&: 

d12  = C.[llFCOS©  4 2 Mp[x3l-COS© -X|PSINq]aX3  j 
f3 3" 

The  fourth  component  of  the  g column  rector  is: 

^ "*■  "^44  ^ . . . ( A . 1 ?d ) 

where 

e ,2  = <e>  I3p  ^ilN  © 

f44=  -G 

A scrutiny  of  tie  preceding  term  by  term  developments 
will  disclose  that  the  option  of  rigli-fln  dynamics  is 
obtained  by  dropping  terms  involving  the  q3  and  qA  gener- 
alized coordinates  and  their  corresponding  coefficients. 

The  " A " matrix  reduces  to  a 2x2  matrix  in  the  rigid-fin 
situation  as  an  illustrative  example. 
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APPENDIX  B 


Fln-3ynamlcs  Computer  Software: 

The  computation  of  fin-opening  dynamics  as  a time 
history  was  numerically  determined  using  the  Continuous 
System  Modeling  Program  (CSMP)  as  documented  by  Specuhart 
and  Creen,  Ref.  3.  The  numerical  treatment  for  solving  sets 
of  second  order  differential  equations  is  similar  to  the 
analog  computer  approach  in  that  the  first  step  is  the 
determination  of  the  highest  order  der i va t 1 ve( s ) followed  by 
succesive  integrations.  Each  integration  step  allows  the 
introduction  of  appropriate  initial  conditions. 

The  equations  of  motion  have  been  represented  concep- 
tually in  matrix  form,  eqn.  4,  for  ease  of  visualization. 
It  will  be  noted  that  the  equations  were  lnertially  coupled 
as  developed  by  the  Lagrange  equation  formulation  technique 
which  is  detailed  in  Appendix  A.  An  analog  computer  solu- 
tion approach  could  accommo' ' te  the  inertial  coupling,  but 
vould  have  difficulties  in  i ■ esenting  the  nonlinear  terms. 
A digital  approach  can  nanile  the  nonlinear  terms,  but  does 
require  that  the  equations  appear  as  being  uncoupled  in  an 
inertial  sense.  The  Inertial  uncoupling,  as  mentioned  in 
the  main  section,  was  achieved  by  premultiplying  the  equa- 
tions of  motion  by  the  Inverse  of  tne  inertial  coefficient 
matrix;  i.e.,  . The  existence  of  the  inverse  can  be 
argued  on  the  basis  of  physical  considerations.  Once  the 
second  time  derivatives  of  the  generalized  coordinates  have 
teen  determined  in  an  uncoupled  form,  it  is  a simple  pro- 
cedure to  obtain  lower  order  derivatives  by  using  a fourtn 
order  Runge-Kutta  integration  scheme. 


The 
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being  user  oriented  to  any  person  with  FORTRAN  programming 
skills.  The  technique  emphasized  simplified  input  data  and 
output  result  statements  respectively.  In  addition,  program 
control  statements  were  almost  a direct  description  of  the 
mathematical  equations  or  physical  variables  of  the  problem. 
The  authors  of  CShP  contend  that  tne  technique  allows  the 
user  to  concentrate  on  the  details  of  the  physical  system 
rather  than  the  usual  concerns  of  numerical  analysis  and 
programming.  The  author  of  this  report  is  in  accord  with 
the  contention.  It  was  a pleasant  surprise  to  find  that  a 
small  investment  of  time  in  learning  the  method  produced 
good  dividends  in  the  form  of  results. 


Figure  B.l  shows  a program  flow  chart  Id  order  to 
illustrate  the  procedural  simplicity.  The  CShP  technique 
Identified  the  solution  flo'  in  phases.  For  the  problem  at 
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Fig.  B"! • Program  Flow  Chart 


hand,  only  t ti 3 "initial"  ani  " dynamic”  pnases  were  required. 
As  one  would  suspect,  the  initial  phase  alloyed  the  entry  to 
the  protlen  of  constants  and  parameters.  Tne  initial  condi- 
tion of  fin-opening  angle  rate  was  normally  selected  as  a 
problem  parameter  while  the  initial  missile  roll  rate  vas 
considered  as  a problem  constant.  There  is  a program  option 
to  repeat  the  calculations  vith  changes  being  made  in  the 
constants.  This  latter  option  aided  considerably  to  program 
flexibility,  ani  was  used  tc  advantage  after  the  program  was 
debugged . • 

The  dynamic  phase  is  recognizable  as  the  portion  of  the 
program  that  is  involved  in  the  computation  of. .the  time  his- 
tories.. There  are  two  options  availible,  a "SORT"  and  a 
"NOSORl"  condition.  The  "NDSDRT"  option  allows  the  use  of 
FORTRAN  program  logic  as  entered  by  tne  programmer  with  the 
computational  steps  being  as  specified.  This  section  can 
draw  upon  subroutines  for  program  support  as  in  more  conven- 
tional software.  For  this  problem,  the  "NOSORT"  section 
involved  an  interaction  of  all  the  current  variables  at  any 
given  instant,  which  is  recognized  as  an  inherent  trait  of  a 
nonlinear  problem.  The  fourth  order  Runge-Kutta  integration 
vas  performed  in  the  "SORT"  section.  An  inherent  advantage 
of  this  section  was  the  main  feature  of  the  C5  MP  * namely, 
the  ability  to  identify  internally  in  the  program  the 
sequence  of  computation  steps  so  that  a full  dynamic  inter- 
action of  the  problem  can  be  obtained.  The  programmer  does 
not  need  to  be  concerned  as  to  the  sequence  of  program 
statements  in  this  section?  however,  ltttvill  be  noted  in  the 
program  listing  that  a "people-oriented"  sequence  was  used 
in  order  to  make  the  program  easier  to  read. 

Dutnut  formatting  was  easily  accomplished  by  identi- 
fying a "TITLE"  plus  a listing  of  variables  of  interest  by  a 
".PRINT’’,  statement.  Problem  timing  was  established  by  the 
"TIMER"  statement  which  specified  the  maximum  (or  finish) 
time.  the  problem  time  step  increments  ani  the  output  time 
printing  step  intervals.  There  were  several  numerical 
methods  available  for  the  integration  scheme,  with  the 
choice  defending  u^on  the.  particular  problem.  This  was 
specified  by  the  ’METHOD"  statement.  The  "RANOE"  statement 
provided  a listing  of  the  extremes,  both  minimum  and  maxi- 
mum, for  selects!  variables.  Finally,  the  option  was  insti- 
tuted of  considering  the  problem  as  being  completed  (FINISH) 
when  the  fin-opening  angle  rate  changed  sign,  which  corres- 
ponded in  a physical  sense  to  tne  plastic  collapse  of  the 
crushable  stop  prior  to  fin's  blowing  back  about  its  pivot. 

A program  listing  is  provided  in  Appendix  C in  order  to 
illustrate  the  readability  (and  hence  simplicity)  of  tne 
program.  The  FORTRAN  statements  in  the  main  program  were 
similar  to  usual  techniques  except  that  comment  cards  needei 
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an  asterisk,  * , in  the  first  colu-nn  in  contrast  to  the  "c" 
found  in  the  subroutine  section.  In  addition,  the  DI'IENSIOV 
card  required  a virgule , / , in  column  l and  continuation  to 
a following  card  was  indicated  by  at  the  end  of  tne 
line.  These  distinctions  are  clearly  explained  in  Ref.  zz, 
hence  the  first-tine  user  cf  CSMP  should  not  experience 
great  difficulty. 

Some  difficulty  was  encountered  in  setting  up  the  tine 
interval  -for  the  problem  in  order  tc  maintain  a balance 
between  the  fact  that  increasing  the  tine  step  -night  degrade 
the  tine  history  solution  accuracy  and  -nake  the  ^rnblen 
unstable  and  the  other  constraint  that  too  snail  a tine  step 
interval  could  lead  to  computer  tine  wastage.  A problen 
tine  interval  of  2.5  milliseconds  was  used  at  first,  hut  it 
lei  to  problen  instabilities  following  tne  introduction  of 
the  crushable-stop  logic  into  the  program  during  the  devel- 
opment period.  A time  step  interval  of  2.25  niiliseconds 
worked  quite  well  and  avoided  the  above  instability  without 
requiring  excessive  computational  time.  In  its  final  form, 
the  program  could  provide  a complete  fin-opening  dynamics 
solution  in  about  six  (5)  seconds  per  condition,  a time 
vhich  was  not  considered  excessive. 

Output  results  will  be  found  as  a typical  listing  in 
Appendix  C.  The  case  shown  is  for  an  initial  missile  roll 
rate,  P2(DYi  =),  of  52.33  rad/sec  and  an  initial  fin-opening 
angle  rate,  T2(DY2  =) , of  12.2  rad/sec.  It  will  be  clear 
from  the  listing  that  the  fin  opening  process  was  stopped  by 
the  crushabl  e-s  t Oi>  model  at  approximately  34  milliseconds 
for  the  specified  condition.  Computational  results  are 
discussed  in  the  main  body  of  this  report. 
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APPENDIX  C 


Fin- Dynamics  Program  Listing: 
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ELASTIC  — FIN  OPENING  DYNAMICS 


APPENDIX  D 


Fin- Dynamics  Output  Sample: 
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APPENDIX  E 


Estimation  of  Fln-Openlng  Rate: 

The  release  of  the  fins  by  the  fin  latch  when  the  ER2P 
leaves  the  gun  barrel  is  dependent  upon  a relative  velocity 
taking  place  between  the  latch  and  the  missile  fins.  The 
relative  motion  occurs  from  several  physical  mechanisms 
including .that  the  ER5P  has  a different  deceleration  manner 
than  the  fin  latch  block,  possibly  due  to  the  latch  block 
being  totally  immersed  in  the  gun  muzzle  blast  field  immed- 
iately upon  exiting  from  the  barrel. 

The  latch  is  physically  arranged  to  keep  the  fins 
locked  in  a folded  position  until  the  latch  moves  forward 
along  the  symmetry  aiis  of  the  missile.  A sketch  of  the 
configuration  is  shown  on  Fig.  E-l  for  sake  of  clarity.  The 
forward  motion  of  the  latch  is  associated  with  a linear 
momentum,  and  it  is  the  conversion  of  a portion  of  this 
linear  momentum  into  fin  angular  momentum  which  may  be 
considered  as  a source  of  the  fin's  initial  opening  angle 
rate. 


The  considerations  which  ensue  are  applications  of 
momentum  principles  from  enginnering  mechanics  such  as  des- 
cribed in  Soldstein's  text,  Ref.  1.  The  fin  latch  has  a 
"jump"  or  discontinuity  in  velocity  (at  time  - 2)  due  to  the 
momentum  transfer,  and  the  velocity  ''jump”  in  combination 
with  the  latch's  mass  provides  a measure  of  the  momentum 
transfer.  The  time  period  involved  in  toe  momentum  transfer 
is  extremely  short  and  the  whole  process  may  betconsidered 
somewhat  in  the  manner  of  the  application  of  a "Dirac  delta" 
type  of  forcing  function.  The  velocity  "jump"  is  defined 
by : 


av  = [V (0+)  - V (3-)]  ...(  E.l  ) 

For  purposes  of  the  ensuing  analysis,  a unit  value  of  vel- 
ocity "jump"  corresponding  to  one  foot  per  second  will  be 
assumed  In  order  to  provide  an  estimate  of  fin-opening  rate 
sensitivity.  Pertinent  dimensions  for  tns  development  are 
shown  on  Fig.  E-l  for  sake  of  clarity,  even  though  some  are 
duplicated  in  Table  I. 
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V V 

^-Fiw  Lpcvch 

(UWC  DWQ.  SK  4 5033 57) 


^-Fin  (tor6) 

(NV4C  DW6.  SVC  45039  2*0} 


C.G.  LOCATION 
X,  =-0.774  in. 
X5—  5.225  in. 

L/VTCH  COHTA,CT 

X,=~  1.770  in. 

X*  = M.00  IN. 


Fig.  E -l : Sketch  of  Fin  Latch 
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It  vlll  be  assumed  that  the  impulse  from  the  fin  latch 
acts  along  the  X3axis  in  the  forward  direction  (-X3  ).  The 
radial  impulse  contribution  is  assumed  to  be  reacting  inter- 
nal to  the  latch  since  all  six  fins  share  equally  in  the 
radial  component.  Therefore,  the  impulse  at  the  latch  con- 
tact point  for  one  fir.  may  be  estimated  as: 

<^3i=  A I ■ (1/6)  HUAV  ..♦(  E.2  ) 

'•  = -2.003314  lb-sec 

where 

» - weight  of  fin  latch  = 2.64  lb 

Ml  = mass  of  fin  latch  = 2.221657  lb-sec*-in' 

AV  = velocity  "jump"  = 12  in-sec-1 


The  radial  contribution  to  the  impulse  transfer  is  based 
upon  the  true  impulse  acting  normal  to  Lne  fin  trailing  edge 
cutout  “ici  a 12  degree  slope  as  sketched  on  Fig.  E-l . 


= <7’3l/(  tan  12°)  = 2.2197?  lb-sec 


i t 


.3  ) 


In  solving  for  the  initial  opening  rate,  three  momentum 

- * i e t V.  rt  Q ^ n ^ ) 4*  A c I n ^ D a r \ Tk  ( » » 1 1 C.O 

LUU  5CA  * O H y ll  L‘4UU  V IVUd  WU  j w <»•---  r 

reaction  at  the  fin  ^ivot  is  involved  In  the  process. 
Conservation  of  linear  momentum  (X  , direction): 

5^p.  - M £ cssN'  ©O)  * & •••(  E-4  ) 

Conservation  of  linear  momentum  ( X 3 di recti  on } : 

^ + M Jl  sin  'H'  GO)  = 2 • • • ( F,.5  ) 

Conservation  of  angular  momentum  (about  fin  c.g.): 


U 0(3) 


3i  u 


' w. 


- xa) 

- x,  ) 


♦ ) 


. ..(  E, 
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Id  the  above  conservation  expressions,  It  is  reoognl2ed 
that : 


V - fin  weight  = 2.93  lbs 

M = mass  of  a fin  = 2.227591  lb-sece-in-1 

I 2 = fin  mass  moment  of  inertia  component  about 
the  c.g.  = 2.09358  lb-in-sec2 

A = (X,2  + x/)'/2  = 5.2£2  in. 

sin  ^ = 2 . 1465 

cos  ^ = 2 .9392 

Substituting  numerical  values  into  eqns.  E.4  through  E.5 
yields : 

2.23966  ©(2)  = -0.01679  lb-sec 

- 2. id 537  ©(2)  = +0.3  2331  lb-sec 

5 . 225  + 2.774  2 . 29358  0(  2 ) « +0.12521  lb-sec 

Finally,  the  numerical  solution  for  initial  fin  opening 
angle  rate  becomes: 

6(3)  = 0.659  rad-sec'  ...(  E.7  ) 

which  is  a result  of  an  impulse  input  from  the  fin  latch  due 
to  a unit  velocity  "jump"  of  one  foot  per  second. 
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